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So far it has been shown that the quantum dynamics cannot be described as a classical Markov 
process unless the number of classical states is uncountably infinite. In this Letter, we present a 
stochastic model with time-correlated noise that exactly reproduces any unitary evolution of a qubit 
and requires just four classical states. The invasive updating of only one bit during a measurement 
accounts for the quantum violation of the Leggett-Garg inequalities. Unlike in a pilot-wave theory, 
the stochastic forces governing the jumps among the four states do not depend on the quantum 
state, but only on the unitary evolution. This model is used to derive a local hidden variable model, 
augmented by one bit of classical communication, for simulating entangled Bell states. 



It is a well-established fact that the quantum dynam- 
ics among a finite set of mutually exclusive alternatives 
(like the up and down states of a 1/2-spin) cannot be re- 
duced to a Markov process among a finite set of classical 
states. Indeed a time-homogeneous Markov process on a 
finite space always relaxes into a stationary probability 
distribution, as a consequence of the Perron-Frobenius 
theorem [l|. The impossibility of exactly simulating a 
quantum system through a Markov process on a finite 
classical space is also a consequence of a theorem proved 
in Ref . ■ There we showed that the dimension of the 
classical space cannot be smaller than the quantum state 
manifold dimension under the hypothesis of Markov dy- 
namics, that is, the classical space must be uncountably 
infinite. Pilot wave theories provide an example of this 
overflow of classical resources. Indeed in their framework 
the dynamics among the set of alternatives explicitly de- 
pends on the quantum state, which actually turns to be 
part of the classical description. Pilot wave theories of 
finite-dimensional quantum systems were considered for 
example in Refs. [1,13] • While they reintroduce a classical 
realistic picture of the quantum world, they are unavoid- 
ably characterized by a feature that is absent in prequan- 
tum physics, namely, the invasiveness of measurements. 
Thus, measurements do not provide a mere updating of 
knowledge about the actual state, but intrinsically in- 
troduce a perturbation on the system. In every known 
classical model of quantum dynamics this perturbation 
demands an invasive updating of an uncountably infinite 
amount of information. For example, in the case of a 
qubit, two continuous real variables need to be updated. 

The Leggett-Garg inequalities provide a useful test 
for deciding if a set of data can be explained by a 
measurement-noninvasive classical theory (BJ. Indeed 
they are violated by quantum mechanics. These in- 
equalities are analogous to the Bell inequalities for the 
Einstein-Podolsky-Rosen (EPR) experiment |6| with the 
local measurements being replaced by two consecutive 
measurements on a two-state quantum oscillator. This 
analogy suggests an interesting question. On the one 
hand the violation of the Leggett-Garg inequalities de- 
mands the invasiveness of measurements in any classical 



theory of two-state quantum oscillators. On the other 
hand, the violation of the Bell inequality implies that 
a classical simulation of Bell correlations requires some 
communication between the parties. It is known that a 
finite amount of classical communication, namely one bit, 
can actually account for this violation [7| ■ Thus it is nat- 
ural to wonder if it is possible to simulate the quantum 
dynamics by a classical model that needs the invasive 
updating of a finite amount of information. 

Generalizing the case discussed in Ref. Q , in this Let- 
ter we consider the scenario where two consecutive mea- 
surements of the same observable are performed, at times 
to and ii > io, on a qubit undergoing a generic unitary 
evolution. The outcome of each measurement is one of 
two orthogonal states, denoted by | ± 1). They can be, 
for example, the left and right states of a double-well 
system [5|. We show that this scenario can be classi- 
cally simulated by a two-bit classical model with time- 
correlated noise. The key ingredients of our model are 
time correlation of the noise and minimal measurement 
invasiveness. The first ingredient is required for circum- 
venting the constraint on the classical space dimension 
proved in Ref. 0], as this constraint does not hold for 
non-Markov processes. The second ingredient accounts 
for the invasiveness of measurements by demanding that 
only one bit is invasively updated. Unlike in a pilot-wave 
theory, the noise governing the dynamics is independent 
of the measurement times to an d t\ and depends only on 
the unitary evolution. In particular, the noise is indepen- 
dent of the quantum state. This feature and the finite 
number of classical states distinguish our result from the 
previous ones, such as the class of hidden variable theo- 
ries considered in Ref. p|, where the stochastic matrices 
are supposed to depend on the quantum state. At first 
glance, the finiteness of the classical space seems para- 
doxical since the quantum state space is infinite. Indeed, 
the quantum information is not encoded into the clas- 
sical state of a single execution, but into the statistical 
behavior of many executions. 

Like in a pilot-wave theory, in the presented model 
the system is supposed to be, at any instant, definitely 
in one of the two orthogonal states | ± 1). More pre- 
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cisely, the model contains a two-value discrete index, 
s = ±1, that determines the outcome of a measurement 
on the basis {| — 1), |1)}. If at some time the index s is 
equal to some value sq and a measurement is performed, 
then the outcome is \sq). Furthermore, the measurement 
does not change the subsequent value of the index s. 
The role of s is similar to the role played by the po- 
sition variables in the de Broglie-Bohm mechanics. To 
derive the stochastic model, we first introduce a simple 
measurement-noninvasive model that captures some fea- 
tures of a qubit. The qubit is described by just the bit 
s that is kicked by a time-correlated noise depending on 
the unitary evolution (first ingredient). Since this model 
is measurement-noninvasive, it satisfies the Leggett-Garg 
inequalities and it is not equivalent to quantum mechan- 
ics. We then modify the model by introducing another bit 
that is invasively updated by the measurements (second 
ingredient). The dynamics of s is ruled by both the noise 
and the additional bit. It is shown that this minimal im- 
provement is sufficient for exactly reproducing the quan- 
tum transition between two consecutive measurements. 

It is useful to represent the quantum states by Bloch 
vectors. The unitary evolution is described by a rotation 
on the Bloch sphere. We denote by R(t a , t b ) the rotation 
operator along any time interval [ta,tb] and by vectors 
±ri the states | ± 1). The quantum probability of having 
|si) at time t\, given \sq) at time to, is 

•Pq(si.*i|sQ)*o) = ^[1 + sis n ■ v], (1) 
v = R(t ,tx)n. (2) 

Let us introduce the measurement-noninvasive model. 
First, we define the noise and the rule governing the 
jumps of s. Then, we derive the transition probability of 
s along a time interval. The noise variable is a unit vector 
x(t) that is a function of time. We denote by p[x(t)] the 
marginal probability distribution of x(t) at time t and by 
p{x) the probability distribution of the function x [i.e., 
p(x) is a functional]. The statistical distribution of the 
noise, namely p(x), is defined by the equations 

p[x(t a )} = (4ir)-\ (3) 
x(t b ) = R(t a ,t b )x(t a ) (4) 

for any t a and t b . The first equation gives the marginal 
probability distribution of x(t a ) at time t a . The second 
equation establishes a deterministic relation between the 
value of the noise variable at different times. The pro- 
cedure for generating each realization of the noise is as 
follows. First, we generate a random vector x(t a ) at some 
time t a according to Eq. ([3j. Then, we determine x(t) 
at any time by using Eq. The noise is clearly time- 
correlated, that is, the correlation function (xi(t a )xj(t b )) 
is not equal to zero for t a =/= t b . Given the noise function 
x(t), we need a rule for the dynamics of s(t). We employ 
the simplest deterministic rule by assuming that the in- 
dex s(t) undergoes a jump whenever n ■ x(t) changes the 
sign, that is, whenever x(t) crosses the geodesic of the 



Bloch sphere lying on a plane orthogonal to n. Let us 
summarize the noninvasive model. 

Model 1. Let the unit vector x(t) be a stochastic func- 
tion of time, whose statistical property is given by Eqs. (0 
In each Monte Carlo execution, the two-value index 
s(t) undergoes a jump whenever x(t)-n changes the sign. 

Thus, the Monte Carlo procedure for generating the 
value of s(ti) = s\ at time t\ given s(to) = so a t time to 
is as follows. A noise function x(t) is generated according 
to Eqs. J3J SI). If the signs of x(t ) ■ n and x{t\) ■ n are 
equal (even number of jumps), then s\ is set equal to 
so, otherwise s\ is set equal to — sq. The procedure is 
repeated for each realization. 

Notice that all we need to know about the noise is 
the real function £ : t — ¥ n ■ x(t). In fact, we could 
just suppose that £ is the noise and regard x as an in- 
termediary tool for mapping each unitary evolution to a 
statistical distribution of £, namely, for generating the 
map R — > p(£) from the function R : t a , t b — > R(t a , t b ) to 
the probability distribution (which is a functional) of the 
noise £. As with v and v(t), if not differently indicated, 
we denote by £ the noise function and by £(i) its value at 
time t. Thus, -F(£) is meant as a functional of the noise at 
every time and F[t;(t)] as a function of the noise at time 
t. Notice that the process £ is not Markovian, that is, 
the marginal probability p(£ c |£&>£a) of having £ c = £(t c ) 
at time t c given = £,(t b ) and £ a = £(t a ) at previous 
times t b and t a is not equal, in general, to the marginal 
probability p(£ c |£h)- 

Let us denote by P^(s±;ti\so;to) the transition proba- 
bility from so at time to to s± at time t±, given a noise 
realization £. The rules defining model 1 imply that 
P^(si; ti\so; to) = 6*[siSo« • x(ti)n ■ x(to)], being the 
Heaviside function. Using Eqs. ^ and ((H), we have that 

f<;(si;ii|so;£o) = 0[sis o n- x(ti)v - x{ti)\, (5) 

From Eq. (J5J and the statistical distribution of x(t\) de- 
fined by Eq. ([3]), we find, by the marginalization over 
x{t\), that the probability of the transition so — > si is 

P(si,ti\s ,t ) = ± $ d 2 x0[siso(n- x){v ■ x)] = 

1 i t - -a (6) 

1 — ± arccos(siso« • v), 

Thus, the model does not exactly reproduce the quantum 
probability given by Eq. ([T]). For example, in the case of 
Rabi oscillation between states | ± 1), the quantum prob- 
ability is a cosine squared function of u>{t\ — to), u> be- 
ing the Rabi frequency. Conversely, the model presented 
here gives a triangle function. In particular, for a small 
evolution time the classical probability scales as t\ — to, 
whereas the quantum probability scales as (ii — to) 2 - 
While the model is not exact, it has the nice property 
of generating an oscillatory dynamics, which a Markov 
process on a finite set of states fails to give. This prop- 
erty is granted by the time correlation of the noise. 

Just as a local model satisfies the Bell inequalities, this 
measurement-noninvasive model satisfies the Leggett- 
Garg inequalities, which are violated by quantum sys- 
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FIG. 1: Schematic representation of the 4-state model. Dur- 
ing a unitary evolution the index s undergoes jumps whenever 
£ r changes sign, whereas r remains constant. For example, if 
a measurement is made at time to and r is set equal to —1, 
the jumps occur at the filled circles. 

tems. To simulate exactly the quantum transition be- 
tween two measurements we need to introduce some ad- 
ditional variable that is invasively updated by the first 
measurement. We now present an exact model that uses 
just one additional bit, denoted by a discrete index, r(t), 
taking the values ±1. The model of a qubit is as follows. 
Model 2. Let the unit vectors x\{t) and X-\(t) be two 
stochastic junctions of time. They are statistically inde- 
pendent and the statistical property of each function is 
given by Eqs. JS^- The qubit is described by two indices 
s{t) = ±1 and r(t) — ±1, which are functions of time t. 
The index s(t) is directly measurable at any time and is 
not modified by a measurement, whereas r(t) is invasively 
updated. If a measurement of s(to) = so is performed at 
time to, then the index r(to) is set equal to 

r = sign {[xi (t ) • nf - [£_i(t ) • n} 2 } , (7) 

The index r(t) remains constant after the measurement, 
whereas s(t) undergoes a jump whenever x ro (t) -ft changes 
the sign. A second measurement reveals the value of 
s(ti) = si at time t\ > to. 

Like in the previous model, all we need to know about 
the noise are the functions £±i : t — > H ■ x±i(t). In the 
following we will denote by £ the pair of noise functions 
A schematic representation of the model for a par- 
ticular realization of the noise is given in Fig. [TJ Notice in 
figure that a measurement at time to sets r = — 1, since 
(,-i(to) > £i(*o)) in accordance with Eq. (|7|). 

For each noise realization £, the probability of having 
outcome s\ at time t\, given outcome so at time to, is 
equal to P^(si;ti|s ;to) = 0[sis o ft ■ x rn (ti) ft ■ x ro (t )] 
where ro is given by Eq. (J7]). Thus, using Eq. ([2]) and 
Q, we have that 

Pjt(si;ii|s ;io) = 9[sis n-x ro {h) v ■ x ro (ti)}. (8) 

Similarly, Eq. ((7]) can be written, through Eq. ^ and 
©, as 

r - sign {[(fx (h) ■ v\ 2 - [£_i(ti) • v} 2 } , (9) 

Like in the previous model, the probability 
P{si; t±\so; to) for the transition from sq a t time to 



to si at time t\ is obtained by averaging over the noise 
realizations. Thus, from Eqs. Q3I8I9P we have that 

P(sr,ti\s ;t ) = (i^E^i/^i^-i . 
9(siSo ft ■ x r v ■ x r )9{r[(xi ■ v) 2 — • v) 2 }}. 

Noting that the two terms in the sum over r give the 
same contribution, it is not difficult to show that 

P(si\ti\so',to) — J d 2 xd(siSo ft ■ xv ■ x)I(v ■ x) (11) 

with I{n) = d 2 y9[r 1 2 - (v ■ y) 2 ] = ±\n\. This gives 

the equation 

P(si;ti\s ;t ) = J d 2 x9(s n- x)p ks (x;siv), (12) 

where pks{x] w) = ■ x9(w ■ x) is the probability distri- 
bution associated with quantum state w in the Kochen- 
Specker (KS) model Q. The integral in Eq. ((TJ) is well- 
known [8| and gives the quantum probability of having 
the state soft, given the state siv and viceversa. Thus, 
we have proved that 

P(si;ti|s ;to) = PQ(si;ti\s ;t ), 

that is, the stochastic model exactly reproduces the quan- 
tum transition between two measurements. Unlike in a 
pilot-wave theory, the noise £ does not depend on mea- 
surement times t\ and ti. Thus, the noise value £(t) at 
any time t is independent of the quantum state at that 
time. Indeed the information on the quantum state is 
encoded in the correlation between (s,r) and £. 

There is a close relation between this model and some 
results in quantum communication. Let us consider the 
following EPR scenario [g. Two spatially separate par- 
ties, Alice and Bob, each receive one of two maximally 
entangled qubits. Alice performs a local projective mea- 
surement on the single-qubit basis (vi,— v*i), while Bob 
on the basis (vo, —Vq). The two- value indices so and si 
arc defined so that Alice and Bob's outcoming states are 
the Bloch vectors —sivi and so^o, respectively. With 
a suitable choice of the reference frame on the Bloch 
spheres, the joint probability distribution of s and si is 
P e {so,si\vo,vi) = j[l + sosi^o • According to Bell's 
theorem, a local hidden variable model cannot reproduce 
this probability distribution and some post-measurement 
communication between the parties has to be exchanged. 
How much communication is required? In Ref. @ it was 
shown that a finite amount of communication, namely 
8 bits, is sufficient for reproducing the Bell correlations. 
This result was improved in Ref. [7] , where it was shown 
that an exact simulation demands a communication of 
just one bit. An alternative model with minimal commu- 
nication was derived in Ref. [l(| from the Kochen-Specker 
model [§[. The common setting of a one-way classical 
protocol for simulating entanglement is as follows. Bob 
and Alice share a random variable X. Given the mea- 
surement (vo,— vo), Bob generates the outcome sq and 
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an additional discrete index m according to a probability 
distribution that depends on v and the shared variable 
X. Then he sends m to Alice. Alice generates the out- 
come si of the measurement (Hi, —V\) with a probability 
that depends on both v\ , X and the communicated index 
m. A stochastic model of quantum dynamics on a finite 
classical space, such as that introduced in this paper, can 
be easily converted into a model of entanglement, where 
the stochastic noise and the classical state play the role 
of X and the communicated information, respectively. 

Let us show that the stochastic model of qubit de- 
rived here can be converted into the model of entangle- 
ment reported in Ref. [lOj. Suppose that at time to the 
qubit is in the mixture |(|1)(1| + | — 1|) and a pro- 
jective measurement is performed on the basis (ft, — ft). 
The probability distribution of outcome so is p(sq) = h. 
Then the qubit undergoes two consecutive unitary evo- 
lutions along the time intervals [to,t] and [i, ii]. At time 
t\ another measurement on the basis (ft, —ft) is made 
and it gives outcome si. Each measurement at time 
ti, with i = 1,2, is actually equivalent to a measure- 
ment on the basis (vi,—Vi) at the same time t, where 
Vi = R(ti,t)ft. Thus, we have that the joint probability 
of so and si is formally equal to the probability distri- 
bution P e (so, Si\vo,vi) for two entangled qubits in the 
EPR scenario. Indeed, both the marginal distributions 
and the correlations of sq and si are identically repro- 
duced. From our stochastic model we find that, given 
the vectors Xi(t) = y*i, the joint probability distribution 
of (s, r) at time t and outcome sq at time to is 

PB(s ,s,r\y;v) = 
i0[( So «-^)(^-y r )]0{r[(«.y 1 ) 2 -(«.y_ 1 ) 2 ]}, [ ' 

where y = (y±, y~i). Similarly, we have that, given y and 
(s, r), the probability of the outcome s± at time t\ is 

PA(si\s,r,y;vi) = 6[(s n ■ y T ){s x v x ■ y r )\. (14) 

Finally, the joint probability of s an d S\ is 

P e (so,si\v ,vi) = 
J2 s ,r I d A yP A (si\s, r, y; vi)P B (s ,s, r\y, vo)p(y). 

These three equations give a model of entanglement 
where y is the shared noise and (s,r) the communi- 
cated bits. Notice that ft is just a free parameter of 
the model. It can be eliminated by the transformation 



s — > s sign(ri-y r ). Furthermore, the marginal probability 
distribution of s after the transformation is uniform and 
independent from v. Thus, s can be included in the set of 
shared variables and, indeed, incorporated in y. In this 
way we obtain a local hidden variable model of entangle- 
ment, augmented by one bit of communication (namely, 
r) , given by the conditional probabilities 

Pb(sq, r\y; v Q ) = 6(s v Q ■ y r )0{r[(v Q ■ y x ) 2 - (v Q • j/-i) 2 ]}, 

Pa(si \r, y;v\) = 0(siih- y r ). 
This is the model derived in Ref. [Hi, set in a slightly 
different form. The process can be reverted and one can 
obtain the stochastic model directly from the model of 
entanglement. 

In conclusion, we have presented a stochastic model 
with time-correlated noise that exactly reproduces any 
unitary evolution of a qubit by using just 4 classical 
states. The time correlation of the noise allowed us to 
overcome the constraint of the theorem proved in Ref. 
on the dimensionality of the classical space. A general- 
ization to higher-dimensional quantum systems can have 
some interesting implications. First, it would automati- 
cally give a local hidden variable model of entanglement, 
augmented by a finite amount of one-way communica- 
tion, for n Bell states. Apart from some approximate 
protocols reported in Ref. [1(| , such a model is at present 
missing. Second, this generalization can suggest more 
efficient methods for simulating the dynamics of high- 
dimensional quantum systems. At first glance, this ap- 
proach does not seem to provide any computational ben- 
efit. Indeed, the evaluation of each noise realization re- 
quires one to solve the Schrodinger equation, thus it is 
not less complicated than directly solving the dynamics 
of the quantum state. However, we have seen that only 
partial information about the noise (in our model, the 
real functions £j) is actually involved in the dynamics of 
the discrete indices. Thus, one could envisage a com- 
putational strategy for computing, exactly or with some 
approximation, this partial information without passing 
through the Schrodinger equation. 
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